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1. Introduction 

Spherically symmetric self-similar perfect fluid spacetimes have attracted considerable 
attention during the last couple of decades (see e.g., [jl], || for references). This is not 
surprising since they lead to ordinary differential equations while still providing an 
arena for a number of interesting physical phenomena: shock waves Q ; self-similar 
perturbations of the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) solution, 
which are of relevance for the growth of primordial black holes and evolution of voids 
(see, e.g., [[jj || ||); and violation of cosmic censorship 0, are some examples. These 
models also have broader implications. For example, they have turned out to be of 
essential importance when it comes to understanding spherically symmetric black hole 
formation j|, ||] . 

Self-similar spherically symmetric perfect-fluid models exhibit several preferred 
geometric structures. There are preferred directions associated with the 4-velocity of 
the fluid and the homothetic Killing vector (with self-similar spacetimes here we mean 
spacetimes with a homothetic Killing vector. For a discussion about self-similarity, 
see ]l0|). Another preferred structure is the area of the spherical symmetry surfaces. 
These three structures have led to three approaches, each using a specific coordinate 
gauge. The first is 'the comoving perfect-fluid approach' initiated by Cahill and Taub 
[B|. The second is 'the homothetic approach' used by Bogoyavlensky and co-workers 
[I. 11. [l^]. The third is 'the Schwarzschild approach' used by, for example, Ori 



and Piran (see, e.g., 0). Each gauge suggests more or less natural variables and 
each gauge has its individual physical interpretational advantages. Thus they are all 
complementary. However, it is possible to express many features gauge invariantly. 
This makes it possible to algebraically relate the variables one has used in the various 



approaches (this will be done in Appendix B|) . It is worth mentioning that apart from 



the above three approaches there are also some other approaches, useful for studying 
the spacetime structure 'near' singularities and for investigating global spacetimc 
features, based on synchronous coordinate systems and null coordinates ||13| . 

We will use the homothetic approach as our starting point. The homothetic 
approach has the advantage that it reveals a considerable structural similarity 
between the field equations for the self-similar spherically symmetric models and 
the hypersurface homogeneous models (e.g. spatially homogeneous (SH) models). 
There is a considerable literature about how one can treat the field equations of the 
hypersurface homogeneous models. Thus the homothetic approach makes it possible 
to transfer ideas from the hypersurface homogeneous arena to the present situation. 
However, the homothetic approach, where one chooses a coordinate along the orbits of 
the homothetic Killing vector (see, e.g., |l4|, [ll| 16 ), has a disadvantage. In general, 



the symmetry surface changes causality. If one uses a homothetic diagonal gauge 
for the spherically symmetric case, then the spacetime must be covered with several 
coordinate patches; one for the region where the homothetic Killing vector is spacelike 
and one where it is timelike. Then one has to join the regions where the homothetic 
Killing vector is null. 

This is the first paper in a series with the aim of obtaining a global picture 
of the space of solutions for self-similar spherically symmetric perfect-fluid models. 
Here we will consider the case where the homothetic Killing vector is spacelike; the 
so-called spatially self-similar (SSS) spherically symmetric models. These models 
are characterized by a four-dimensional homothetic symmetry group acting multiply 
transitively on three-dimensional spatial surfaces. In a subsequent paper p7| we will 
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deal with the timelike self-similar (TSS) spherically symmetric perfect-fluid models, 
i.e. the models where the homothetic Killing vector is timelike. Most of the physically 
interesting phenomena, e.g. shock waves, are actually associated with the TSS region, 
and consequently it is not surprising that most previous work has focused on this 
region. However, to obtain a global understanding one also needs to consider the 
spatial region. In a third paper we will address the specific problem of matching 
solutions between the two regions as well as focusing more on physical properties of 
the space of solutions. Therefore the mathematical nature of the present paper is 
justified as the results are needed in the subsequent studies. However, one physically 
interesting property pertaining to the SSS region, density perturbations of the flat 
FLRW solution, is treated in the present paper. 

By using the homothetic approach and ideas from hypersurface homogeneous 
models we will reformulate Einstein's field equations so that the dynamics takes place 
in a reduced and compact phase space; one system for the spacelike case and one for the 
timelike case. Moreover, it will turn out that in our formulation, all equilibrium points 
(critical points, singular points, fixed points) are hyperbolic (here we also include the 
case of a line of equilibrium points where there is one zero eigenvalue corresponding to 
the line). This is in stark contrast to earlier treatments where noncompact variables 
have been used and where certain parts of phase space have been 'crushed'. To deal 
with these problems, one has made numerous variable changes and special treatments 
of various parts of phase space (see, e.g. This has resulted in an incomplete, or 
even misleading, picture of the global solution structure. 

The line element for SSS spherically symmetric models, written in diagonal form 
where one of the coordinates is adapted to the homothetic symmetry, takes the form 
1 

ds 2 = e 2x ds 2 = e 2x [-dt 2 + D\{t)dx 2 + D 2 2 {t) (d6 2 + sin 2 {9) dip 2 )] .(1) 

We will consider perfect fluid models. The energy momentum tensor, T ab , is thus 
given by 

f ab = jlu a u b + p{u a u b + gab), (2) 

where fx is the energy density, p is the pressure and u a the 4-velocity of the fluid. We 
will assume 

P = (7-1)A. (3) 

as an equation of state, where the parameter 7 takes values in the interval 1 < 7 < 2 
which include radiation (7 = |). Consequently, we have excluded dust (7=1) and 
stiff fluids (7 = 2). These models behave quite differently compared to those in 
the above interval, and thus need special treatment. The dust solutions are known 
explicitly (these models are just special cases of the general spherically symmetric dust 
solutions, which are all known explicitly, see e.g., p8|). 

The outline of the paper is the following. In section 2 Einstein's field equations 
are rewritten in terms of a dimensionlcss set of variables in order to obtain a maximal 
reduction of the coupled system of ordinary differential equations. The variables are 
chosen so that they take values in a compact phase space. In section 3 the equations 
are subsequently analysed by means of methods from the theory of dynamical systems. 
In section 4 a numerical investigation is undertaken and global dynamical features are 
considered. We end with a discussion in section 5. Appendix A describes the properties 
of the fluid congruence, the condition for the spacetimes to belong to Petrov type 
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and also some expressions relevant for the mass function and density perturbations. 
Appendix B gives the relation between various coordinates and the relation between 
the present variables and other variables used in the literature. 

2. The dynamical system 



Following JTjj] we introduce 
D X = B 1 ~ 1 



o /3 -2/3+ 



Do = Bo 



o /3 +/3 + 



(4) 



3/3' 



bo 



3/3+, 



where the dot denotes d/dt. The quantities 9 and cr+ describe the kinematical 
properties of the normal congruence of the symmetry surfaces in the (M,ds 2 ) 
spacetime that is conformally related to the physical spacetime (M,ds 2 ) with the 
homothetic factor (see, e.g., 9 is the expansion while <r+ describes the shear. 

The orthonormal frame components, in (M, ds 2 ), of the fluid velocity arc conveniently 



parametrized by (1, v, 0, 0)/vT 



where v is just the 3-velocity with respect to 



the symmetry surfaces. Einstein's field equations, G ao — T a b, and the conservation 



equations, T' 



ab 



0, lead to a set of equations presented in |19[. However, for 



mathematical reasons we will now go over to another set of variables. The motivation 
lies in a simplification of the constraint while keeping the 'canonical quadratic' nature 
of the defining equation for [i n (see the equations below). The new variables are given 

by 



1 



(•5) 



-^(29- 



1 

3 



= -={9 + 2a + ). 



This leads to: 



Evolution equations 



1 
1 

V3 



1 



9 2 + a% + 9a+ - 3Bf 
a\ + 29a+ 



qd2 , 3( 7 -l)(l-^ 2 ) 

l + ( 7 -l)« 2 Mr ' 



3(2- 7 ) + (3 7 - 2)v 2 



l + ( 7 -l)i; 2 



B 1 =-^a + B 1 , 

B2 = -^=(9 + a+)B 2 , 
1-v 2 



V3 7 (l- (7- 1)« 2 ) 
+v / 3 [(7 - 1)(3 7 - 2)v 2 - (2 - 7)] Bi} 



{ 7 [2(7-l)£ + 7 a + ] v 



(6) 



Constraint equation 



7«Mn"^ [l + ( 7 -l)f 2 ] a + B 1= 0. 



(7) 
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Defining equation for fx n 



Mn = I (O 2 'SB 2 + 3£?f ) . (8) 



Auxiliary equation 

20+{l-v 2 )a+ + 2V3vB 1 ] . (9) 



2" 



\/3 [1 + (7 - 1) v 

The quantity /i„ is the energy density of the fluid measured by an observer with the 
4- velocity orthogonal to the symmetry surfaces in (M, ds 2 ), see | [L6| . 

As the final step we will introduce a new set of variables where the scale invariance 
is used to decouple one equation. In spatially homogeneous cosmologies one often 
uses the expansion for this purpose. However, a very desirable property is to obtain a 
reduced phase space which is compact. Unfortunately, using, for example, 9 or 9 does 
not achieve this in the present case. Instead we will use the quantity 9 2 + 3B 2 since 
this is a 'dominant quantity' (see equation (||) and note that we assume a nonnegative 
energy density). We define a new variable Y, 



Y = ^9* + 3Bl B 2 = ^J{Y^-9^)/i. (10) 

Unfortunately, this will lead to a more mathematical treatment than if one uses a more 
'physical' quantity instead of Y . However, there seems to be no such quantity which 
leads to a compact phase space. Instead one is forced to algebraically relate interesting 
physical objects, e.g. fluid kinematic quantities and mass functions to the variables 



used in this paper. These relations are given in Appendix A. We also introduce the 
Y"-normalized bounded dimensionless variables 

Oo = |, Q+ = ^, Cx = 4^. (ID 

Note that C% is positive. The density \i n is replaced by the density parameter Qy, 
which is defined by 

An introduction of a dimensionless time variable f 

dt = VSY^df , (13) 
leads to a decoupling of the ^-equation 



Y' = F Y Y, F y 



,(14) 



where a prime denotes d/dr. The remaining coupled evolution equations can now be 
written in dimensionless form: 
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1 (2 - 7) + (37 - 2)v 2 
2 



C[ =2Ci 



1 + (7- l)i; 2 
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(15) 



2 1 + ( 7 -l)w 2 
{ 7 [2(7-l)go + 7Q+] w 



7 [1 — (7 — l)v 2 j 
[(7-l)(37-2)« 2 -(2- 7 )] 

Constraint equation 

G = jvn Y - 2(1 + (7 - 1) v^Q+d = 0. 



Defining equation for fly 
tt Y = 1 - 



(16) 
(17) 



Auxiliary equation 



7 



^y-^ [2g + (1 - w 2 )g+ + 2vCi] + 2F y | .(18) 



l + (7- 

The reduced equations and the equation for Y are invariant under the transformation 

(f, Qo, Q+, C u v) (-f, -Qo, -Q+,Cx, -v). (19) 

The line element can be obtained when Qo, C\ and Y have been found through the 
relations 



Dx = V3(YCi)-\ D 2 = V3[Y 2 (l-Ql) 



-1/2 



* = V^/y -(20) 

Because of the above discrete symmetry, there are, in general, £wo orbits in the reduced 
phase space, related by (|l~9|), corresponding to a single line element. There is also a 
class of orbits that are invariant under the above discrete symmetry and in this case 
each orbit corresponds to a single line clement. 



3. Dynamical systems analysis of the reduced phase space 

The reduced phase space is determined by (Qo, Q+, C\, v), related by the constraint 
G = 0, given by equation (|l6|). The boundary of the physical phase space is of 
essential importance in understanding the dynamics of the interior phase space and 
hence we will include it in our analysis and thus obtain a compactificd phase space. 
We will not solve the constraint globally since this leads to problems. Instead we will 
follow |l^, |2(J and solve it locally around the equilibrium points (i.e. we will solve 
the linearized constraint for different variables at different equilibrium points). This 
formulation will enable us to achieve a good understanding of the global structure of 
the reduced phase space. 
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3.1. Equilibrium points and local analysis 

Here we will present the equilibrium points together with Sly , which will indicate if a 
point is a vacuum point or not, the gradient of the constraint and what variable we 
have locally solved for, and finally the eigenvalues. There arc numerous equilibrium 
points. However, they are often closely related. We will use the following notation: 
the equilibrium points are denoted as 

s g n(Q )Kernei:-g +) . (21) 
When there is no risk for confusion we will suppress sgn(Qo), sgn(Q+) or sgn(v). 

3.1.1. Equilibrium points with zero tilt, v = 0. 
The equilibrium points ±F: 

Qo = ±1, Q+ = -^, Ci = 0; fV = | 

VG = (0, 0, Qo, y ), (Ci eliminated); (22) 

Ai =^(3 7 -2)Q , A 2 = ^(2- 7 )Qo, A 3 = ±(3 7 - 2)Q„. 

These points are saddle points in the full phase space but there is a two-dimensional 
separatrix surface, spanned by the eigenvectors corresponding to the eigenvalues Ai 
and A3 above, the orbits of which lie in the interior of phase space. 

The equilibrium points ±K±: 

Qo = ± 1, Q+ = Qo, C x =0; fl Y = 0; 

VG = (0, 0, -2Q , 0), (Gi eliminated); (23) 
Ax =2Q , A 2 =3(2- 7 )Qo, A 3 = (3 7 - 2)Q . 

The point (-K°_) is a local source (sink). Therefore there is a two-parameter 

set of orbits starting (ending) at the point +K+ (-K°_). 

The equilibrium points ±K^.: 

Qo = ± 1, Q+ = -Qo, Gi - 0; Oy = 0; 

VG = (0, 0, 2Q , 0), (Gi eliminated); (24) 
Ax =2Q , A 2 = -A 3 = (2- 7 )Q . 

These two points are saddle points on the boundary and there are no asymptotes into 
the physical phase space. 

3.1.2. Equilibrium points with intermediate tilt, < v 2 < 1. 
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The equilibrium points ±M : 



Qo = ± 1, Q+ = 0, Cx = 1, 

-7(7- 1) + 



VG 



= 0; 



= 



(7- 



Ai = 



2[1 
-2Q , 
(1 + Q v) [(2 



A 2 = 



l)( 7 2 (7-l) + (2-7)(3 7 -2)) 
(3 7 -2)( 7 -l) 

l)v 2 ] ,-2^v, 0) , (Ci eliminated); 
2Q (7-1)W 



Go; 



(25) 



(3 7 -2) [l-( 7 -l)„2] 
7) + (3 7 - 2)Q v] 



where 



W = (3 7 - 2) + 4(3 7 - 4)vQ + (3 7 2 - 4)v 2 - (3 7 - 2) (7 - l)v 4 .(26) 

These points correspond to the Minkowski spacetime. The constraint |vj < 1 implies 
that these two points are physical only for 7 > |. The point + M (-M) enters the 
physical part of phase space though the point + M + (_M~) (discussed below) at 
7 = |. Note the bifurcation in the eigenvalue for 7 = | in equation (|2^). As 7 
increases, +M (-M) moves along the line + M + — + AI~ (_A/~ — _A/ + ). From the 
expressions above, it follows that A2 < and A3 > for | < 7 < 2. Therefore +M 
(_M) is a saddle point in the full phase space. The positive (negative) eigenvalue A3 
is associated with the vacuum submanifold so there is a two-dimensional separatrix 
surface ending (starting) at +M (-M), the orbits of which lie in the interior phase 
space. 



3.1.3. Equilibrium points with extreme tilt, v 2 



1. 



The equilibrium points ±K± : 

Qo = ± 1, Q+ = Qo, Ci =0; n Y = 0; 

VG = (0, -2 7 sgn(w), -2 7 Qo, 0), {C x eliminated); 

2(3 7 -2)Q 



Ai — 2Qo 7 A2 — 4Qoj A3 



'7 



(27) 



Note that the sign of v affects the linearized constraint but not the stability of the 
points. These points are saddle points on the boundary and there are no asymptotes 
into the interior phase space. 



The equilibrium points ±M^: 

Qo = ± 1, Q+ = 0, Gi = 1, v = Qo; 



n Y = 0; 



VG = (0, -2 7 , -2 7 Qo, 0), (Gi eliminated); 
Ai = — 2Qo, A2 = — 4Qo, A3 = 



(28) 



2-7 

These points correspond to the Minkowski spacetime. Note that there is a bifurcation 
f. Forl< 7 <! 



at 7 = |. For 1 < 7 < I the points ±M ± are saddle points on the boundary and there 
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are no asymptotes entering the interior of the phase space. At 7 = | they coincide 
with the points ±M respectively. For | < 7 < 2 the point + M + (_M~) is a sink 
(source) and there is a two-parameter set of orbits ending (starting) at the point. 

The equilibrium points ±rt^ : 

Qo = ±1, Q+ = -Q (l-C 1 ), v = -Q ; 
n Y = -2QoQ + (l + Q Q + ); 

VG = (0, -2 7 , 2Q07, -2(2 - i)QoQ+{l + QoQ+)); (29) 

(Gi eliminated); 
Ai = A 2 = -2Q (l + 2QoQ+), A 3 = 0. 

The constraint fly > implies that Q+ is restricted to the interval — 1 < Q+ < 
(0 < Q+ < 1) when Qo = 1 (Qo — —1)- There is only one non-zero eigenvalue 
correponding to two different eigenvectors. Equilibrium points on the line +H~ (-W + ) 
are sources when — 1 < Q+ < — | (0 < Q + < \) and sinks when — \ < Q + < 
(|<Q+<1). 

The equilibrium points ±KZ: 

Qo = ± 1, Q+ = -Qo, C x =0; Q Y = 0; 

VG = (0, 27Sgn(u), 2 7 Q , 0), (Gi eliminated); (30) 

Ai =A 2 = 2, A 3 = 0. 

These two points are the endpoints of the physical part of the lines ±H T , 
respectively. They describe the boundary asymptotes which play an important role 
in understanding the global dynamics of phase space. This motivates their special 
treatment. The points have no asymptotes entering the interior phase space. 

The equilibrium points ±M^ : 

Qo = ± 1, Q+ = 0, C 1 = 1, v = -Qo; V Y - 0; 

VG = (0, -2 7 , 2 7 Q , 0), (Gi eliminated); (31) 

Ai =A 2 = -2Q , A 3 = 0. 

These points correspond to the Minkowski spacetime, and are the other endpoints 
(compared to ±K0) of the lines ±H T . These points have no asymptotes entering the 
interior phase space. 

3. 2. Invariant submanifolds on the boundary of the SSS phase space 

The boundary is described by a number of invariant subsets: Qo = ±1, &y — 
1 — Q+ — G\ = 0, Gi = 0, v — ±1. The submanifolds with Qo = ±1 are described by 
the same equations as the reduced equations for the plane symmetric type il models. 
The equations corresponding to the submanifolds with v = ±1 are identical to the 
reduced field equations for SSS spherically symmetric models with a neutrino fluid as 
the source. The submanifolds associated with fly = are vacuum submanifolds with 
a test fluid. The constraint ( fL6| ) implies that Q+Gl = when fly = 0. If Q + = 0, one 
obtains the same equations as one has for SSS spherically symmetric vacuum models 
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Table 1. The various boundary submanifolds. 



Boundary 


Restriction 




ill 


Qo = ±1 




V 


n Y = o, 


. = o, d = 1 


N± 


V = ±1 




KSV± 


Oy = 0, Q^ 


. = ±1, Ci = 


KS 


Ci = 0, v = 






with a test perfect fluid, while C\ = leads to the same equations as those describing 
spatially homogeneous Kantowski-Sachs (KS) vacuum models with a test perfect fluid. 
Apart from the above vacuum submanifolds, C\ = also leads to a submanifold with 
v = 0. This submanifold is described by the same equations as the orthogonal SH 
KS perfect-fluid models. The submanifolds are collected in Table 1 where we also 
introduce designations. 

We will now describe the dynamical features of the individual boundary 
submanifolds. 



The ±\I submanifolds. The dynamical structure of the plane symmetric type +il 
submanifold is given in figure 1(a) (7 > |) and (b) (7 < f ). That of „il is easily 
obtained by applying the discrete symmetry ([l9|). 

The V submanifold. The dynamical structure of the V submanifold is given in figure 
1(c) (7 > I) and (d) (7 < |). Note that Qo is an increasing monotonic function. 

The submanifolds. The dynamical structure of the N~ submanifold is given in 
figure 1(e). That of N + is easily obtained by applying the discrete symmetry (|l9|). 
Note that Q+ < (Q+ > 0) in the N~ (N+) case. The N ± submanifolds are solvable. 
The orbits are characterized by the following integral: 

Q+Ci 



[Qo -Q+- sgn(w)Ci 



= constant. (32) 



The KSV± submanifolds. The dynamical structure of the KSV- submanifold is given 
in figure 1 (f ) . That of KS V+ is easily obtained by applying the discrete symmetry ( |l9| ) . 
The variable v satisfies the inequality v < (v > 0) in the KSV_ (KSV+) submanifold. 
The KSV± submanifolds are solvable. The integral describing the various orbits is 
given by 



[1 - sgn(Q + )g ] 7 



constant. (33) 



It is worth pointing out that Qo is an increasing (decreasing) monotonic function 
and that v is a decreasing (increasing) monotonic function for the KSV_ (KSV+) 
submanifold. 
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The KS submanifold. Since the equations for this submanifold are the same as those 
for the SH KS models we automatically, as a bonus, get a dynamical description of 
these models. The SH KS models are of considerable interest and have been discussed 
in the literature before. Collins formulated the field equations as a dynamical system 
and made a qualitative analysis [plf . However, he did not obtain a compact phase 
space. Uggla and von Zur Muhlcn discussed the locally rotationally symmetric SH 
Bianchi type IX models The KS models appear as a boundary submanifold in 

this context, but this was unfortunately not pointed out in |^]. Neither was a phase 
portrait given. The dynamical structure of the KS submanifold is given in figure 1(g). 
Note that Qo is a monotonic function taking the solutions from big bang to big crunch. 
There exists a single time-symmetric orbit invariant under the discrete symmetry ( |l9| ) . 
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V 



-1 



I 




Figure 2. This picture shows how the boundary of the reduced phase space of 
the spatially self-similar spherically symme tric models is constructed out of the 
invariant submanifolds discussed in section 3.2. 



There is also a single orbit starting (ending) at the point +F (-F). The + F point (_F 
point) corresponds to the flat FLRW solution in the SH KS geometry, and consequently 
these models have an asymptotically flat FLRW behaviour. 

Orbits in the interior phase space coming close to the KS submanifold can be 
approximated with KS orbits. The corresponding self-similar models in the interior 
phase space can thus be approximately described by a SH model multiplied by the 
nomothetic factor, i.e. they are approximately conformally SH. 

4. Global dynamical behavior 

By appropriately 'gluing together' the boundary submanifolds of the previous section, 
we obtain the reduced phase space shown in figure 2, of which the SSS models are 
orbits in the interior. By indicating the stability of the equilibrium points on the 
boundary we get figure 3 for j < | and figure 4 for 7 > |. The orbits on the 
boundary corresponding to eigenvector directions of the points ±F and ±M are also 
shown. The two-dimensional separatrix surfaces entering the interior phase space from 
these two points will be discussed in section 4.3. 

4-1- Monotonic functions 

Monotonic functions are important tools for understanding the global dynamics. For 
example, the existence of a monotonic function in a region of phase space implies that 
there cannot be any limit cycles in that region. 
The function 




(34) 



whose derivative is given by 



M' 



(2- 7 )(3 7 -2)(l- W 2 )C 1 



71/. 



(35) 
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is monotone in the regions v < and v > 0. There are no interior invariant 
submanifolds with v = 0. Thus all attractors (in the present case, equilibrium points) 
lie on the boundary. 



4-2. Classification of orbits 

It is natural to divide the orbits in the SSS phase space into three categories depending 
on the extendibility properties into the timelike self-similar (TSS) region. In order to 
continue an orbit into the TSS region it must start or end at an equilibrium point with 
extreme tilt (\v\ = 1). The classification presented here is based on whether orbits 
start and/or end at ±H T . It is also possible to extend orbits into the TSS region 
through ±M , but from the structure of the TSS part of the phase space |17 it is 
clear that these orbits are not as physically interesting as the orbits starting or ending 
at j-TP. The three categories are: 

0: Orbits that are purely SSS, i.e., orbits that do not start or end at ±H T . 

I: Orbits that may be continued into the TSS region 'once', e.g. orbits that start 
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Figure 5. Separatrix submanifolds in the interior of the reduced phase space, 
a) The +F sub-manifold. The dotted orbit is the FLRW-orbit. b) The -M 
submanifold. c) The S submanifold. The dotted orbit is the TOVKMZ-orbit. 
Dashed orbits lie in the boundary submanifolds. 



(end) at an equilibrium point not belonging to ±TP while they end (start) at ±TC^ . 

II: Orbits that may be continued into the TSS region 'twice', e.g. orbits that 
start and end at ±H T . 

4-3. Invariant submanifolds in the interior of the SSS phase space 

The ±F submanifolds. The separatrix surface spanned by the outgoing (ingoing) 
eigenvectors of the +F (-F) point will be referred to as the + F (-F) submanifold 
(figure 5(a)). The + F submanifold contains two different types of orbits. Near the 
+ il submanifold, the orbits belong to category I and end in a band on the +'H~ line. 
The separatrix boundary in +i I defines the upper boundary of the band. The lower 
boundary is for Q + — — |, where the stability of the +H~ line changes to a source. 
The orbits in the +F submanifold that are close to the Kantowski-Sachs submanifold 
belong to category 0. They approach +7Y~, but are repelled and end at -K^_. The 
structure of the -F submanifold is easily obtained by applying the discrete symmetry 
(^9|). Note that for the + F (_F) separatrix, v < (v > 0) in the whole interior. 

In an SH slicing, an additional nomothetic Killing vector leads to solutions being 
represented as equilibrium points |23||. This is the case with the flat FLRW solution. 
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initial singularity 



Figure 6. Schematic picture of the fiat FLRW spacetime with a slicing adapted 
to the homothetic Killing vector, compared to a slicing adapted to the comoving 
time T. The cones are hypersurfaces of constant self-similar variable t. The region 
inside the lightcone is the TSS region. Note that the cones actually have a much 
more complicated structure than in this over-simplified picture. 



However, in the present case we have a slicing associated with the homothetic Killing 
vector, see figure 6. In this formulation, the flat FLRW solution appears as an one- 
dimcnsional orbit. The SSS part of the flat FLRW solution actually corresponds to 
two (equivalent) orbits contained within the 'category I band'. The FLRW orbits can 
be found by imposing the Petrov type condition along with vanishing fluid shear 
and acceleration. The orbits are characterized by 

2 <2° Cl = - . „ . (36) 



',+ 



4+(37-2)w 2 ' 4+(37-2)w 2 ' 

The FLRW orbit in the + F submanifold ends at a point on +H~ with Q + = 
-2/(37 + 2). 

The invariant quantity 2m /R (see Appendix A. 2) computed for the FLRW orbits 
is given by 

4 

(37) 



2m - 1 r, r 2 



(37 - 2) 2 v 2 ' 

Note that there is a one-parameter set of orbits within the +F submanifold lying 
arbitrarily close to the FLRW orbit throughout the entire SSS region. 

Self-similar fluctuations in a flat FLRW universe have been studied by Carr and 
Yahil Q, although they focus their discussion on phenomena pertaining to the TSS 
region. In order to study density perturbations, we follow Ellis and Bruni ^4j. Their 
approach has the advantage of being covariant and gauge independent. The fractional 
density gradient is characterized by a quantity L (see Appendix A. 3). This quantity 
vanishes identically for non-tilted SH models. In the present case we have two such 
examples, namely the flat FLRW solution and the entire KS boundary submanifold. 
Solutions close to the flat FLRW in the + F submanifold can be parametrized by L: 
negative values of L correspond to over-dense solutions (p, > Pflrw), and positive 
values correspond to under-dense solutions. Thus, the flat FLRW solution splits 
the +F submanifold into an under-dense part and an over-dense part. Over-dense 
solutions, starting out from the equilibrium point + F in a direction sufficiently close 
to the KS boundary, behave completely differently compared to the FLRW solution 
and correspond to category orbits. The remaining orbits in the +F submanifold 



Spatially self-similar spherically symmetric perfect-fluid models 



16 



behave roughly as the flat FLRW orbit and belong to category I. Near the equilibrium 
point + F, the fractional density gradient vanishes to linear order. Consequently, 
all solutions in the + F submanifold are near-homogeneous initially, with density 
fluctuations growing significant only at later stages, possibly with the exception of 
the origin r = 0. 

The ±M submanifolds. The equilibrium points ±M exist when 7 > |. In this case 
one has additional separatrix surfaces in phase space. As for the + F submanifold, 
the orbits in the _M submanifold (figure 5(b)) are separated into two bands. Orbits 
starting close to the vacuum submanifold belong to category I and end on +7i~ in 
a continuous band from + M~ down to Q + = — i. The other band contains orbits 
belonging to category 0, which are repelled from +H~ and end at -K°_. 



The Symmetric submanifold. Orbits passing through (Qo — 0, Q+ — 0, C\.v = 0) 
constitute a submanifold containing 'time-symmetric' solutions (figure 5(c)). The 
orbits in the symmetric submanifold (S submanifold) fall into two classes: those 
starting at -7i + and ending at +7Y~, thus belonging to category II, and those 
belonging to category 0, starting at +K+ and ending at -K^_. 

One of the orbits in the S submanifold belonging to category II corresponds to 
the SSS part of a static self-similar solution associated with many names: Tolman 
p5| , Oppenheimer and Volkoff |26| , Klein |27j] and Misner and Zapolsky [^8) , see e.g. 
Henrikscn and Patel |l3|, Collins 29 and Bogoyavlensky Q. We will call this orbit 
the TOVKMZ orbit. It is characterized by 

== 2( 7 -l)Qo = 7Qo r „ s s 

Q+ - 37^' Cl -~(37^2K (38) 

The quantity 2m /R takes the form 

~R -(2- 7 )' + 8( 7 -l)- (39) 



4-4- General structure of interior orbits: 

Here we consider the behavior of orbits, not necessarily belonging to any of the 
submanifolds discussed above. A general orbit must start at one of the sources of 
the reduced phase space, and end at one of the sinks. The analysis is performed 
separately for 7 < | and 7 > |. For 7 < |, the possible sources are the point +K^_, 
and the 'source parts' of the lines +H~ and _7Y + . The only sinks are the point -K^_, 
and the 'sink parts' of the lines +H.~ and _7Y + . Thus there are nine possible types of 
orbits for 7 < |. It seems highly unlikely that orbits starting at +H.~ and ending at 
-7i + exist, especially since the -K^_ sink in general is very strong. Apart from such 
orbits, the existence of general orbits of all categories have been verified numerically. 
When 7 > |, the appearance of the ±M points changes the stability of _Af~ to a 
source, and + M + to a sink. Apart from all the types of orbits occuring for j < |, there 
will be additional types of orbits. However, some of the new possible combinations 
of start- and end-points do not occur. This is due to the ±M separatices 'screening' 
the ±M^ points. Thus, orbits starting at _A/~ and ending at + M + or _7Y + do not 
exist. Likewise, orbits starting at +TL~ and ending at + M + do not appear either. 
Apart from these exclusions, the existence of general orbits of all categories have been 
verified numerically for 7 > | as well. 
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5. Discussion 

In this paper we have used a dynamical systems approach to study the global structure 
of the solution space of SSS spherically symmetric perfect fluid models. By choosing 
variables adapted to the homothetic Killing vector field the phase space of these models 
can be made compact - a very desirable property. As the constraint equation for 
these models is not solved globally, but rather locally at the equilibrium points of the 
system, we avoid 'crushing' parts of phase space. As a result of having a complete 
regular phase space we find that there is a bifurcation at 7 = |, something which has 
been overlooked in the previous literature. Also, a compact and regular phase space 
allows us to follow trajectories throughout their entire evolution. 

It is of interest to observe that the FLRW model appears as a trajectory in phase 
space. This is in stark contrast to Bianchi cosmology where the flat FLRW appears 
as an equilibrium point. This affects the perturbation analysis of this model. In 
Bianchi cosmology one only needs to know the phase space in a neighbourhood of the 
equilibrium point. In the present case the FLRW trajectory spans most of the phase 
space and thus there is a need to have a full phase space description to carry out a 
perturbation analysis. 

The existence of a monotonic function implies that there are no limit points 
in the interior of phase space. This stresses the importance of the boundary. The 
SH Kantowski-Sachs models appear as a boundary submanifold of the reduced phase 
space, and thus we obtain a compact description of these models as well. 

We see that in the phase space corresponding to the SSS region there are models 
which are close to the FLRW orbit throughout their entire evolution in the SSS 
region. The models corresponding to these orbits can be interpreted as cosmological 
models with density perturbations. There are also orbits arbitrarily close to the static 
TOVKMZ orbit which also can be continued into the TSS region. To obtain a full 
spacetime picture of the models corresponding to these perturbed FLRW and static 
orbits, we will have to extend them into the TSS region. The TSS region will be 
studied in a forthcoming paper. 

Appendix A. Interesting physical quantities 

Appendix A.l. Fluid properties and Petrov type conditions 
The fluid expansion is given by 




e 



— X 



V3(l - v 2 f/ 2 



{vv' + (1 - v 2 )(2Q a + Q+ + 3v(7i)} Y. 



(A.l) 



The fluid shear is given by 



a = 




e 



— X 



3(1 - w 2)3/2 



{vv — 9 



2a + } 
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e~ x 



3(1 - W 2)3/2 

The fluid acceleration is given by 

3v 



{vv'-(l-v 2 )(Q a + 2Q + )}Y. (A.2) 



3vr 



0v - 2cr + u + 3Bi 



+ 3V3Bi| 



73(l„ w 2)3/2 

The Weyl scalar is 



3^3(1 - v 2 f/ 2 

{v'~ {l-v^iQ+v-C^Y. (A.3) 



C 2 = C abcd C abcd = le- ix (3<7 + - 3B 2 2 + 0a+ - 2a\) 



4 -4x 
9 e 



2 



!) 



Vz(6 + 2&+\ - 3B| - 6a+ - 2a 2 
4 e- 4x [Q' + 2Q' + + F Y (Q + 2Q+) 

+Ql - Q+ (Qo + 2Q+) - 1] 2 Y 4 . (A.4) 

Note that the magnetic part of the Weyl tensor is identically zero for all models. 
The spacetime is of Petrov type D if C 2 7^ and Petrov type if C 2 = 0. 

Appendix A.2. The mass function 

For general spherically symmetric spacetimes the total mass-energy m between the 
centre distribution and some 2-spacc of symmetry is defined as, see e.g., Misner and 
Sharp (30) , Hernandez and Misner Q and Cahill and McVittie Q , 

2m _ a6 OR 8R 1 + 2Q Q Q+ + Q\ - C 2 

R 9 8x a dx b l-Ql ' 1 j 

where the parameter R is the invariant radius of the symmetry surfaces. By definition, 
the matter is outside the gravitational radius 2m whenever 2m/ R < 1 and inside the 
gravitational radius when 2m / R > 1. The latter case is associated with the existence 
of an apparent horizon located at 2m/ R = 1. 

Appendix A.3. Density perturbations 

A covariant approach to density perturbations has been given by Ellis and Bruni pij . 
The fractional density gradient is defined as 

^) = pr^dbii. (A.6) 

f* J a 

In the reduced phase space variables, this becomes 

A a =Ye- x —^(-v, 1,0,0), (A.7) 
1 — v z 



where 



_ 2C 1 [1 + (7 - IK] + iv [2Qo + (1 + v 2 )Q+] 
v / 3[l-( 7 -l)« 2 ] 



and Y is the decoupled variable (10). 
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Appendix B. Coordinate and variable transformations 

Here we will give transformations to other coordinates and variables which have been 
used in the literature to study self-similar spherically symmetric models. The present 
variables lead to a compact and regular description everywhere in the SSS region. This 
is not the case with previously used variables, as is easily seen in the equations below. 
The simple algebraic relations make it easy to identify the points, or even manifolds, 
where the breakdowns occur. 



Appendix B.l. The comoving (fluid) approach 

The most commonly used approach is the fluid approach, where one uses a coordinate 
system adapted to the fluid velocity. The line element, of the presently considered 
SSS spherically symmetric models, takes the following form in the fluid approach (see, 
e.g., @): 

= _ e *(A) dT 2 + e A(X) dx 2 + Y 2 (\)X 2 dQ. 2 , (B.l) 

where A = X/T and dfl 2 = d8 2 + sm 2 (6)dtp 2 . The fluid velocity, u a , is given by 
e-*/ 2 (l,0,0,0) ; 

The following transformation: 

X = e x - F ®, r = e B - F © +f I (B.2) 

where 

dF 
~dT ' 

leads to the 'homothetic' adapted SSS spherically symmetric line element 

d~s 2 = e 2x ds 2 = e 2x [-N 2 dP + D, 2 dx 2 + D 2 2 dSl 2 ] , (B.4) 

where 

N 2 = e A+*+K-2F (e A _ eW+af)-^ 

D\ = e- 2F (e A - e* +2f ) , D 2 = e~ 2F Y 2 . (B.5) 

The fluid velocity v is given by v 2 = e* +2i ~ A . Inserting this into the above expressions 
we find that 

dF v 2 „ v 2 e A - 2F 



(B.3) 



at 1 — 1 — v 2 

The fluid approach was used by Foglizzo and Henriksen see also Bicknell and 
Henriksen |(J , to write the self-similar spherically symmetric equations in a very simple 
way. Their variables are related to those used in this paper according to 

3^ 2 v 2 n n 3jva+ = 3jvQ+ 

K! !i 4[l + (j-l)v 2 ]B 2 2Ci 



[l + (j-l)v 2 ] (9 + a+) 2 + 3B 2 - 3Bf 



l±F-H 



(1 - V 2 ) fl n 

V3-/V US + a+f + 3B\ - 3B\ 



2Sia+(l - v 2 ) 
3jv (1 + 2QoQ+ + Q 2 +~ C\\ 



(B. 



2CiQ+(l-« 2 ) 

vf-h =v~\ (B.9) 
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Appendix B.2. The Schwarzschild approach 

The usual way to represent a spherically symmetric line element is of the form 

ds 2 = -FdT 2 + GdR 2 + R 2 dVL 2 (B.10) 

where, in the self-similar case, F and G are functions of R/T only. The following 
transformation: 



where the function <f> satisfies 
d(j> _ V3B 2 



T : 



c+<p(t) 



C 2 Y 



(B.ll) 



(B.12) 



dt 9 + a+ V3(Qo + Q+)' 
transforms the line element equation ( |B.10| ) into the homothetic form. The metric 



functions F and G can be written as 



F 



3(Q + Q+) V 2 *Y" 



B\ 3B 2 - (6 + cr 4 



G = 



3B 2 



C 2 [C 2 - (Qo + 



a + ) 2 C 2 ~(Q + Q+) 2 ' 



3B 2 - 

Identifying 2m/ R = 1 — G _1 leads to equation ( |A.5| 
The radial 3-velocity vr of the fluid is given by 

V3vBi +9 + d+ 



(B.13) 
(B.14) 



VR 



vCi + Qo + Q+ 
V3Bi +v(6 + a+) ~ d + v (Qo + Q+) 



(B.15) 



In their treatment of self-similar spherically symmetric collapse, Ori and Piran j/j used 
the following set of variables: 

In - \2 



M 



3B 2 -3B{ 1 + 2Q Q+ + Q+ - C{ 



R 



6Bi 



D = 



(1 - v 2 ) fx Tl 



2 (1 - Ql) 



,(B.16) 



2(l + ( 7 -l)«3)^ 



2jv(l - Q 2 ) 



VSvBx 



cr + vCi + Qo + Q+ 



V3B 2 Vl - v 2 
The variables of Maison M are just 



(l-Qg) (!-«*) 



A 2 



B AI 



Pm 



3B 2 - (0 + a + f _ C 2 - (Qo + Q+) 



3B 2 



1-08 



3B 2 



<?+) 2 {Qo + Q+f 



$.17) 
(B.18) 

(B.19) 
(B.20) 



4irR 2 fl 



(1 - v 2 ) Mr 



2B 2 [1 + ( 7 ~ 1)« 2 ] 

(1 -U 2 ) fiy 

2 (1 - Qg) |l + (7-lH' 



VM = VR. 



(B.21) 
(B.22) 
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